Exchange Interaction and T c in Alkaline-earth-metal-oxide-based DMS without 
Magnetic Impurities: First Principle Pseudo-SIC and Monte Carlo Calculation 
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The prospects of half-metallic ferromagnetism being induced by the incorporation of C atoms 
into alkaline-earth-metal-oxides are investigated by the first principle calculation. The origin of 
the ferromagnetism is discussed through the calculation of the electronic structure and exchange 
coupling constant by using the pseudo-potential-like self-interaction-corrected local spin density 
method. The Curie temperature (T c ) is also predicted by employing the Monte Carlo simulation. It 
is shown that by taking the electron self-interaction into account, the half- metallic ferromagnetism 
induced by C in the host materials is more stabilized in comparison with the standard LDA case, 
and the C's 2p electron states in the bandgap become more localized resulting in the predominance 
of the short-ranged exchange interaction. While the ferromagnetism in MgOi-ajC^ is stabilized due 
to the exchange interaction of the lsf-nearest neighbor pairs and might be suppressed by the anti- 
ferromagnetic super- exchange interaction at higher x, the ferromagnetism in CaOi-id, SrOi-^C^, 
and BaOi-iCi is stabilized by both the 1st- and 2nd-nearest neighbor pairs, and T c monotonously 
increases with the C concentration. 
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Besides the attempts to discover dilute magnetic semi- 
conductors (DMSs) by incorporating transition metals 
into various materials to realize ferromagnetic DMSs for 
spintronic devices, recently the ferromagnetism induced 
by nonmagnetic impurities has also attracted great atten- 
tion among both theoretical and experimental scientists. 
A novel class of magnetic materials can be formed by in- 
corporating nonmagnetic impurities or by lattice defects. 
With regard to this issue, the magnetism induced by the 
cation vacancy in MgO [1| , Ca vacancy in CaO > Hf 
vacancy in Hf0 2 with T c exceeding 500 kHHHT^ 

and 

by hydrogen in graphite [H, [^] and in carbon nanotubes 
[l(j | etc. has been reported. It has also been predicted 
that magnetism can be induced by nonmagnetic impu- 
rities such as C and N by substituting the O atoms in 
various oxides. For example, the ferromagnetism caused 
by N and C in alkaline-earth-metal-oxides [H, [ll| and 
the half-metallic ferromagnetism induced by N in quartz- 
Si0 2 0. 

It is suggested that the magnetism might arise in the 
host materials if impurities have a finite local magnetic 
moment that interacts with each other to form a mag- 
netic moment net. In the materials in which the mag- 
netism is induced by incorporating nonmagnetic impuri- 
ties, the substitutional ions may have a nonzero magnetic 
moment and the 2p-electrons of these ions, rather than 
the 3<i-electrons, play an essential role in inducing the 
magnetism in the host materials. They form an impurity 
band in the deep bandgap, and ferromagnetism can be 
induced if the Fermi level lies in these impurity bands. 
Furthermore, recent research also shows the role of the 
2p-like impurity band formed by C and N in the stabi- 
lization of the ferromagnetism in Gai-^Mn^As [l3j and 
Im_ x Mn x N M- 
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In this Letter, we will discuss the origin of the ferro- 
magnetism that is induced by substituting C for O in four 
alkaline-earth- metal-oxides — MgO, CaO, BaO, and SrO 
— through the calculation of not only the electronic struc- 
ture but also the exchange coupling constant and Curie 
temperature (T c ). It should be noted that the prediction 
of ferromagnetism and T c in ref. 3 and ref. 11 is based 
on the electronic structure calculation using the density 
functional theory within the local density approximation 
(LDA) and the mean field approximation (MFA) which 
cannot be successful for describing many materials, es- 
pecially for strongly correlated systems. The LDA often 
overestimates the hybridization between electron states 
due to the underestimation of the bandgap energies of 
semiconductors. Moreover, the calculation of T c within 
the MFA often predicts T c with an excessively high value 
even if the substitutional concentration is lower than the 
percolation limit fl5l ]. Therefore, in order to describe and 
predict the properties of these systems, we need a more 
accurate method. There have been several attempts to 
improve the LDA to overcome these inaccuracies. One 
of the most popular methods is LDA+U which improves 
the LDA by taking the strong correlation effects into ac- 
count through the screened Coulomb parameter U [H[. 
The other popular approaches are the self- interaction cor- 
rection (SIC) [ItJ and the recent pseudo-SIC [HI meth- 
ods. SIC methods differ from LDA+U in that LDA+U 
uses the additional parameter U, while the SIC methods 
do not require any additional parameter at all. Con- 
sequently, in order to investigate the origin of ferromag- 
netism in alkaline-earth-metal-oxide-based DMSs, we ap- 
ply the pseudo-SIC method based on the improvement 
of the MACHIKANEYAMA2000 package coded by Akai 
[l9| (LDA+SIC) in the calculation of the electronic struc- 
tures, and then we use the formula by Liechtenstein et 
aL[20| to calculate the exchange interaction Jy between 
two impurities at the ith- and jth-sites in the ferromag- 
netic coherent potential approximation medium. Finally, 



2 



i 

i 

E 1 
i. 

■ 

* r 
2 , 





1 




1 C Jto'DOB 






\ 


/ ? 



■ Ire**) 

FIG. 1: Total DOS (solid line) and C's 2p - PDOS (dashed 
line) for CaO0.9C0.1- The left-hand figure (a) shows the stan- 
dard LDA calculation and the other (b) shows the pseudo-SIC 
results. The upper plane corresponds to the majority spin and 
lower plane to the minority spin. 



we employ the Monte Carlo simulation to estimate T c . 

The substitution of O with C in alkaline-earth-metal- 
oxides is treated randomly. For convenience, the lattice 
constants of alkaline-earth-metal-oxide-based DMSs are 
fixed to the values of undoped crystals [2l[ (a = 4.123 A, 
4.909 A, 5.160 A, and 5.520 A for MgO, CaO, BaO, and 
SrO, respectively) and no distortion in lattice structure 
is assumed. Throughout the electronic structure calcu- 
lations, 624 independent fc-sampling points in the first 
Brillouin zone are used. The potential form is restricted 
to the muffin-tin type, and muffin-tin radii are chosen in 
such a way that the ions at lattice sites can touch each 
other. The relativistic effect is also taken into account 
by using scalar relativistic approximation. 

In order to compare the density of states (DOS) 
calculated within the standard LDA with that within 
LDA+SIC we plot the DOS of the typical case of alkaline- 
earth- metal-oxides, CaOi-xCir, at x = 0.10 in Fig. Q] 
The figure on the left-hand side corresponds to the DOS 
calculated within the standard LDA (Fig.[lja)), and the 
figure on the right-hand side illustrates the DOS within 
LDA+SIC (Fig. QIb)). It can be seen that, taking the 
self-interaction of electrons into account, the LDA+SIC 
calculation gives a wider bandgap than the standard 
LDA. While the position of the minority spin states of 
C's 2p electrons remains unchanged, the majority spin 
states calculated within LDA+SIC shifts about 1.1 eV 
in comparison with the standard LDA one, leading to a 
higher localization of LDA+SIC 2p states than LDA. And 
the local magnetic moment of C increases from 1.242/1 ^ 
(LDA) to 1.482/xs (LDA+SIC). The exchange splitting 
in the case of LDA+SIC increases approximately 2 times 
as compared to that in the case of LDA, resulting in the 
possibility of the suppression of super-exchange interac- 
tion and the enhancement of the ferromagnetic double 
exchange mechanism. 

Figured depicts the DOS of AOi-xC^ (A = Mg, Ca, 
Sr, Ba) calculated within LDA+SIC at x = 0.10. As seen 
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FIG. 2: Total DOS (solid line) and C's 2p - PDOS (dashed 
line) obtained by pseudo-SIC calculation for MgOo.gCo.i (a), 
CaOo.gCo.i (b), SrOo.gCo.i (c), and BaOo.gCo.i (d). 



from Fig. [21 all materials in question are half-metallic. 
However, the exchange mechanism causing the ferromag- 
netism is somewhat different. Figure [Ha) illustrates the 
DOS of MgO0.9C0.1- With the smallest lattice constant, 
MgOo.gCo.i has the largest bandgap energy. The C's 2p 
states are located near the top of the valence band orig- 
inated from anion p-states, resulting in the strong hy- 
bridization of 2p electron wave functions. The majority 
spin states hybridize with O's 2p states, leading to the 
appearance of the impurity band that connects to the top 
of the valence band and causes a narrower majority spin 
bandgap. The minority spins create an impurity band 
in the bandgap, which includes the Fermi level. This 
impurity band is broadened with a half-width of about 
1.5 eV. The exchange splitting in this material approxi- 
mates to 2.14 eV and is the smallest compared with the 
three remaining materials. These impurity bands can 
be broadened more strongly with increasing C concen- 
tration, leading to the antiferromagnetic super-exchange 
interaction being easy to occur and the compensation 
of majority and minority spins. These mean that the 
ferromagnetism might be suppressed at higher C concen- 
trations. In addition, our calculations for higher C con- 
centrations show that the ferromagnetism in MgOi-zCx 
is most stabilized at x « 0.10 (also refer to Fig. [5]). 
Figures H;b)-[2i;d) demonstrate the DOS of CaOo.gCo.i, 
SrOo.gCo.i, and BaOo.gCo.i, respectively. As seen from 
the figures, the bandgap becomes narrower and the local- 
ization of the 2p states becomes stronger with increasing 
the distance between atoms (or lattice constant). The 
majority spin states are located in the valence band and 
cause a small broadening of the band. The Fermi level lies 
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FIG. 3: Local magnetic moment of substitutional C atoms 
in MgOi-^C^, CaOi-xC^, SrOi-^Ci,, and BaOi-^Ca; at 
x = 0.10. 
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in the impurity band induced by minority spins with 1/3 
of this impurity band being occupied by electrons. The 
exchange splitting in these materials is approximately 2.3 
eV, but the C's 2p states in BaOn.gCo.i is the most lo- 
calized. This change in the localization with respect to 
the lattice constant is also indicated by the value of the 
local magnetic moment of C which is shown in Fig. [3] As 
a result, one can expect the ferromagnetism to be stabi- 
lized by the predominant double exchange mechanism in 
these materials even at higher C concentrations. 

The exchange interactions in MgOi-^C^, 

CaOi-^C^, SrOi-^Cz, and BaOi-^C^ are shown 
in Fig. [U As seen from Fig. BJa), the exchange inter- 
action of the lsi-nearest neighbors in MgOi_ x C x is 
the strongest. The interaction strength is considerably 
strong (about 140 meV for 5% of the C concentration) 
and most of the contributions come from the nearest 
neighbor interaction Jq\. Except Jni and J04 (at 10% 
of the C concentration), contributions from other Jij 
are very small and ignorable. Thus, the exchange 
interaction in MgO\^ x C x can be considered as a typical 
case of short ranged interactions. In addition, the 
exchange interaction is considerably suppressed as the 
C concentration increases. At x = 0.20, the exchange 
interaction becomes very small and the magnetism is 
fully suppressed. Interestingly, by increasing the lattice 
constant, the role of the 2nd-nearest neighbors (J02) 
becomes more important. For CaOi-xCa; (Fig. [5b)), 
Jni is approximately two times smaller than that in 
MgOi-^Cz, but the contributions J02 from the 2nd- 
nearest neighbors becomes important. The dominance 
of the contributions J02 of the 2nd-nearest neighbor- 
pairs in SrOi-^Ca; and BaOi-^Ca, is shown in Fig. (He) 
and Fig. Ud), respectively. In these former materials, 
the exchange interactions of the 2nd-nearest neighbors 
are more important than the lsi-nearest neighbors. This 
can be caused by the sufficiently large distances between 



FIG. 4: Exchange coupling constant vs. distance between C 
atoms in units of the lattice constant a in MgOi-zC^ (a), 
CaOi-^Cz (b), SrOi-zCi (c), and BaOi-zC;,; (d) at several 
concentrations of C (as = 0.05, 0.10, 0.15, and 0.20). 



the lsi-nearest neighbors, while the interactions between 
the 2nd— nearest atoms easily accur through the mediate 
atoms such as Sr or Ba which have large ionic radii. 

In short, the role of the nearest neighbors in the ex- 
change interaction changes with the lattice constant. The 
higher the lattice constant, the greater is the contribu- 
tion from the 2nd-nearest neighbors in the induction of 
the magnetism in the host materials. The exchange inter- 
action is very short ranged and the ferromagnetic double 
exchange mechanism is expected to be predominant for 
all materials in question. Except MgOi-^Ca,, the mag- 
netism in CaOi-zCa,, SrOi-^Cx, and BaOi-^Ca; can be 
more stabilized at higher C concentrations. This is also 
consistent with the results of the Monte Carlo simulation 
for the estimation of T c . 

In order to evaluate T c , we perform the Monte Carlo 
simulations. The Metropolis algorithm [22j is applied 
to calculate the thermal average of the magnetization 
M and its powers. Then, the cumulant crossing method 
proposed by Binder [22j is employed and the fourth order 
cumulant U4 (a linear combination of <M 4 >/<M 2 > 2 ) 
is calculated as a function of temperature for different 
cell sizes (14 x 14 x 14, 16 x 16 x 16, and 18 x 18 x 18 
conventional fee cells) to find the universal fix-point at 
T c . We estimate T c for four values of the substitutional 
C concentrations: x = 0.05,0.10,0.15, and 0.20. The 
obtained results are demonstrated in Fig. [5] Some points 
should be clarified here. 

First, as discussed above for Fig.^a) and Fig0|a), the 
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FIG. 5: Curie temperature T c vs. substitutional impurity 
concentration in MgOi-^d (solid line), CaOi-zC;,; (dashed 
line), SrOi-xCx (dashed-dot-dot line), and BaOi-iCi 
(dashed-dot line). The circles, crosses, and triangles denote 
the calculated points. 



ferromagnetism of MgOi-^C^ is the most stabilized at 
x w 0.10 and will be suppressed at higher C concentra- 
tions due to the anti-ferromagnetic super-exchange inter- 
action and the suppression of the majority and minority 
spins. T c of MgOi-^Ca; (solid line) increases with x in 
the range of x from 0.05 to 0.10 and has a peak of 108 K 
at x 0.10, and then it sharply drops as a function 
of — £ at higher x. As shown in Fig. HJa), the impu- 
rity band formed by C's 2p electrons is strongly broad- 
ened. The bandwidth W of the impurity band in the gap 
can be much larger than the effective correlation energy 
U = E(N+1)+E(N-1)-2E(N), and the Stoner's condi- 
tion for the existence of magnetism might be contravened 
[ill [l2| . Correspondingly, the magnetism might be fully 
suppressed at the C concentrations higher than 16% be- 
cause of the strong broadening of the impurity bands of 
both majority and minority spins, resulting in the com- 
pensation of spins; therefore, the material in question has 
a nonmagnetic behavior at higher C concentrations. 

Second, contrary to MgOi-^Ca;, the impurity states in 
CaOi-^Cj;, SrOi-^Cz, and BaOi- x C x are localized and 
the bandwidths of the impurity bands formed by minority 
spins in the gap are sufficiently small to satisfy Stoner's 
condition; hence, the ferromagnetism is stabilized by a 
predominant ferromagnetic double exchange mechanism 
and is more stabilized as x increases. Although the ex- 



change interaction Joi between the 1 si-nearest neighbor 
pairs considerably decreases with the increasing lattice 
constant, owing to the contributions J02 from the 2nd- 
nearest neighbors, T c monotonously increases with x and 
can gain a value higher than room temperature if the C 
concentration is high enough. However, a question that 
arises here is the solubility of C in the materials. 

In conclusion, we have presented the results of the 
study on the origin of ferromagnetism and predicted T c 
of alkaline-earth-metals-oxide-based DMSs without tran- 
sition metal elements. The electronic structures and ex- 
change coupling constants are calculated by applying the 
pseudo-SIC approach. The dominant exchange mecha- 
nism and the role of the nearest neighbors in these DMSs 
are discussed, and T c is also evaluated by employing the 
Monte Carlo simulation. In short, some comments can 
be made as follows: 

(i) All DMSs in question have a half-metallic ferro- 
magnetism. The ferromagnetic double exchange mecha- 
nism is predominant for all materials in question, except 
for MgOi-^Ca; when the C concentrations is higher than 
10%. 

(ii) The exchange interaction in these materials is 
strong but short ranged. While the contributions 
come mostly from the 1 si-nearest neighbor pairs in 
MgOi-^Cj;, the important role played by the 1st- and 
2nd-nearest neighbor pairs in CaOi-^Ca;, StOi^ x C Xi and 
BaOi-zCz is comparable. However, while Joi becomes 
weaker, J02 gets considerably stronger with the lattice 
constant being larger and becomes the dominant con- 
tribution to the stabilization of the ferromagnetism in 
SrOi-^Cz and BaOi-^Ca,. 

(Hi) Correspondingly, T c of CaOi-^Cz, SrOi^Cz, 
and BaOi- x C x increases monotonously with the increas- 
ing C concentration, while T c of MgOi-^Ca; gains the 
highest value at 10% of the C concentration, then it drops 
sharply and tends to reach zero at higher concentrations 
due to the anti-ferromagnetic superexchange interaction 
and the compensation of majority and minority spins. 
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